Asymmetric Reversal in Inhomogeneous Magnetic Heterostructures 
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Asymmetric magnetization reversal is an unusual phenomenon in antiferromagnet / ferromagnet 
(AF/FM) exchange biased bilayers. We investigated this phenomenon in a simple model system 
experimentally and by simulation assuming inhomogeneously distributed interfacial AF moments. 
The results suggest that the observed asymmetry originates from the intrinsic broken symmetry of 
the system, which results in local incomplete domain walls parallel to the interface in reversal to neg- 
ative saturation of the FM. Magneto-optic Kerr effect unambiguously confirms such an asymmetric 
reversal and a depth-dependent FM domain wall in accord with the magnetometry and simulations. 

PACS numbers: 75.25,+z, 75.60.-d, 75.70.-i 



Exchange coupling between a ferromagnet (FM) and 
an antiferromagnet (AF) has been intensely studied due 
to the fundamental interest in inhomogeneous magnetic 
systems and its central role as a magnetic reference in 
various devices. In most magnetic systems, time rever- 
sal symmetry is present and manifested by a symmetric 
magnetization curve relative to the origin. This sym- 
metry also requires that the magnetization reversal from 
positive to negative saturation be identical to the reverse 
process. However, in a FM/AF system, exchange bias 
(EB) develops below the AF Neel temperature T/v pro- 
ducing a shift (Heb) of the hysteresis loop along the 
magnetic field axis [l|. Therefore, with the shift breaking 
the time reversal symmetry, magnetization reversal sym- 
metry is no longer required. In fact, asymmetric reversal 
was observed by polarized neutron reflectometry , pho- 
toemission electron microscopy 3], magnetotransport Q, 
magneto-optical indicator film 5], and magneto-optical 
Kerr effect In some systems the reversal along the 
decreasing branch is dominated by transverse magnetic 
moments which was interpreted as due to coherent mag- 
netic rotation. The absence of transverse moments in 
the increasing branch reversal was interpreted as domain 
wall propagation 0, 0. Different, even opposite scenar- 
ios were also found 0, 0, 0- Despite the well estab- 
lished experimental evidence and proposed theoretical 
models |£§ |l£J, [ll| , the origin of this asymmetry remains 
a controversial and highly debated issue [l2|. This sit- 
uation is further complicated by the lack of knowledge 
of the interface, crystal imperfections, complex FM and 
AF anisotropy energies, and training effect. While these 
factors are important for each individual system, the fun- 
damental connection of the reversal asymmetry to the 
broken symmetry intrinsic in the inhomogeneous system 
is overlooked. 

In this Letter, we have investigated a simple model sys- 
tem using a variety of experimental techniques combined 



with numerical simulations. We establish a critical link 
between this unusual reversal asymmetry with the time 
reversal asymmetry in these systems. Namely, in reversal 
toward the two FM saturated states, the intrinsic asym- 
metry gives rise to different competing mechanisms, thus 
different reversal processes. 

FeF2/(Ni, Py) bilayers were prepared for this study. 
FeF2 is an AF with a Neel temperature T/v = 78 K, and 
a large uniaxial anisotropy K u — 1.35 x 10 4 kJ/m 3 along 
[001] direction, hence can be considered as a model Ising 
system [ijJLLi' with the AF spins frozen along [001] at 
low temperatures [15|. The Ni or Py (NigiFeig) is poly- 
crystalline with a negligibly small crystalline anisotropy, 
except for a sma ll g rowth-induced uniaxial anisotropy 
along FeF2 [001] 16]. This system is thus in close ap- 
proximation with simple theoretical assumptions. 

The bilayer was grown by e-beam evaporation on a sin- 
gle crystal MgF 2 (110) substrate, where FeF 2 (110) grows 
epitaxially untwinned [HJ . Vector vibrating sample 
magnetometry (VSM) of FeF 2 (50 nm) / Ni (21 nm) 
/ Al (7.6 nm) gives simultaneously the in-plane longi- 
tudinal (parallel to the magnetic field) and transverse 
(perpendicular to the magnetic field) magnetic moments 
[lM LlJj. The magnetic field is applied along the FeF 2 
easy axis [001] with a small misalignment that defines the 
sign of the transverse component during reversal 0, ^| . 
Square hysteresis loops are found above T/v along [001] 
|lq . Cooling the sample in a field ji^Hpc = 0.2 T from 
T = 150 K to 15 K yields an EB field ^H EB = -0.1 T 
(Fig. CJi) and virtually no coercivity. Both longitudinal 
and transverse hysteresis loops exhibit a clear asymme- 
try. Starting from positive saturation, the reversal occurs 
with a sharp corner in the longitudinal component and 
an abrupt increase in transverse component to over 75% 
of the saturation magnetization. Then the FM gradually 
approaches negative saturation, evidenced by the long 
tail in both components. A significant non-zero trans- 
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verse component is found even at hqH = -0.5 T. In the 
increasing field sweep, Ni is saturated almost immediately 
after the reversal. The asymmetry of the two FM orienta- 
tions, especially the long tail around negative saturation, 
is key to understanding the asymmetric reversal. 

We modeled the asymmetric reversal process with mi- 
cromagnetic simulations |2£| using a 20 nm thick Ni 
layer with lateral size 500x500 nm 2 , discretized into 
5x5x2 nm 3 cells. The Hamiltonian H of the system is 
given by, 

H=A [(Vm w ) 2 + (Vm 22/ ) 2 + (Vm u ) 2 ]A7 

i£{FM} 

- (K u mL&V + K d m 2 iz AV + H ■ m) 
ie{FM} 

— Jfm/af ^ rhi-Bi, 

i£{Interface} 

where the three summed terms include FM exchange en- 
ergy, FM anisotropy and Zeeman energy, and FM/AF 
interfacial interaction, respectively. The AF is assumed 
to be frozen during the hysteresis cycle, thus its energy 
contribution remains constant and is not considered in 
the Hamiltonian above. H is the magnetic field ap- 
plied along the x axis with 0.5° misalignment similar 
to the experiment. ra$ and AV are the magnetic mo- 
ment and volume of each cell, respectively. The reduced 
moment fhi is defined by fhi = rhi/Ms. We used the 
nearest-neighbor exchange constant A = 3.4 pj/m and 
the saturation magnetization Ms = 494 kA/m for Ni 
IU . The small growth-induced anisotropy of the Ni layer 
is taken into account by a uniaxial anisotropy along the 
x axis with K u = 5 kJ/m 3 obtained from measurements 
along the hard axis above Tjy. The dipolar interaction 
is approximated by a shape anisotropy along the z axis 
(out-of-plane) with K d = -(/i /2)M| = -150 kJ/m 3 , 
which keeps the moments in the sample plane and avoids 
boundary effects. 

The AF is modeled by a monolayer of spatially inho- 
mogeneous frozen moments, Si, exchange coupled to the 
bottom layer of the FM with an adjustable interfacial 
coupling from Jfm/af ~ Jaf — —0.45 meV up to 2Jaf 
[l^ . We introduced AF grains of average size 25x25 nm 2 
to simulate the inhomogeneous interfacial coupling [l^. 
Gi = —aiSf F pj with Sf F — 2, consists of two random 
quntities: denoting the intergrain variation, and pj 
the intragrain variation, varies as 1±0.35 between 
grains, while pj varies as (7±2)% between cells. This 7% 
assumption is based on recent experiments which found 
net frozen AF interfacial moments with about 4% 22] or 
7% coverage that contribute to EB. Crucial parame- 
ters for the simulation include the product of the uncom- 
pensated moment coverage and interfacial coupling, and 
intergrain fluctuation. The former defines the effective 
coupling strength. The latter describes the interfacial in- 
homogeneity modulated over a length scale of the grain 



size (25 nm), comparable with the FM domain wall width 
82 nm. This spatial modulation of Gi leads to an inhomo- 
geneous pinning on the FM, and is essential to explain 
reversal process revealed in the experiment. However, 
the intensity of the modulation is not essential: 20% to 
50% standard deviation in gives similar results. The 
resultant spatial variation of Gi is shown in the inset of 
Figure ^ 
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FIG. 1: (color online) (a) Vector VSM measurement (filled 
symbol) and micromagnetic simulation (empty symbol) of 
FeF 2 (50 nm)/Ni(21 nm) at 15 K after field cooling in a 0.2 T 
field. Both longitudinal (black symbol) and transverse (red 
symbol) components are measured and simulated. The mi- 
cromagnetic simulation was performed assuming the FM in- 
terfacial layer is coupled to spatially inhomogeneous uncom- 
pensated frozen AF spins, whose distribution is shown in the 
inset (500x500 nm 2 ). The different shades of red refer to the 
different magnitude of local uncompensated frozen AF mo- 
ments with white corresponding to zero local density. 

The simulation (Fig. matches well both the longitu- 
dinal and transverse hysteresis loops exhibiting the same 
asymmetry as the experiment 24] . The bottom and side 
view of the FM spin configuration during the hysteresis 
(Figure 0) shows domains evolving both laterally and in 
the depth. In positive saturation, the FM is uniformly 
magnetized throughout the thickness because both the 
applied magnetic field and interfacial coupling favor this 
orientation. As the magnetic field decreases, the reversal 
is initiated from the top of the FM far away from the 
interface while the bottom pinned by the AF remains 
in the positive direction. An incomplete (non-180 ) FM 
domain wall (IDW) is thus formed parallel to the in- 
terface. As the field decreases further, these FM IDWs 
slowly shrink laterally and squeeze close to the interface. 
Even at jhqH = -0.8 T, the FM is not saturated at some 
interface regions. This lateral domain formation is the 
result of the spatially varying Gi . The regions in the FM 
most resistant to reversal are where the strongest local 
interfacial pinning is found. As the field increases, these 
regions become nucleation sites for the development of 
local IDWs both laterally and in the depth. Therefore, 
these local IDWs result from the competition between 
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inhomogeneously distributed interfacial pinning and the 
magnetic field. Due to the unidirectional nature of the 
AF pinning field, it only competes with the Zeeman en- 
ergy in approaching negative saturation, while they both 
stabilize the FM when positively saturated. This simu- 
lation demonstrates that the local development of IDWs 
constitutes the dominant asymmetric reversal mode. Al- 
though similar exchange spring is claimed in hard/soft 
magnetic structures [23, |2£|, it does not lead to asym- 
metric reversal In addition, this incomplete domain 
wall is unusual in EB because the interfacial coupling en- 
ergy is much weaker than that in a conventional exchange 
spring, thus it was never convincingly observed and was 
overlooked in most EB studies. 

When a finite anisotropy of pinned AF moments is in- 
cluded in the simulation, the IDW is pushed into the 
AF forming a hybrid domain wall across the interface, 
but the main features of the reversal process remain un- 
changed. Since the anisotropy of the FM is usually much 
smaller than that of the AF, the FM-side of domain wall 
dominates the reversal. 









11 








J 











FIG. 2: (color online) Images in the first row from left to right 
are the simulated FM spin configurations (500x500 nm 2 ) at 
the FM/AF interface at 0.8, -0.08, -0.36, and -0.8 T, respec- 
tively; the second row shows the simulated FM depth profiles 
(125x20 nm 2 ), the bottom edge referring to the FM/AF in- 
terface) for the same corresponding field cross sectioned at the 
green lines. Red, white and blue corresponds to M x /M s = 1, 
0,-1, respectively, with x being the magnetic field direction. 

The result implies several important features of the lo- 
cal IDW reversal process. First, the FM domain wall 
depth-dependence is crucial for the asymmetric rever- 
sal process. An important signature of this behavior is 
the asymmetric development of transverse magnetic mo- 
ments. This behavior tends to be smeared out by AF 
twinning or polycrystallinity, and/or more complicated 
FM or AF anisotropy energy terms. It is worth not- 
ing that this asymmetry of approaching two saturated 
states may seem different from the asymmetry of two 
field sweep branches observed before in other systems, 
where a sharp corner is found in the decreasing branch 
and a rounded one for the increasing one |j, |6| . How- 
ever, they are essentially the same except for the small 
FM uniaxial anisotropy, thus negligible coercivity in our 
system. If the FM uniaxial anisotropy is increased to K u 
= 50 kJ/m 3 and a 0.5° fanning of the AF pinning moment 



in the sample plane is included, the simulated hysteresis 
loop displays the same asymmetry as observed before to- 
gether with an irreversible transverse loop (Fig. |3J) 
Second, the local nature of the IDW due to the interfacial 
inhomogeneity is crucial in the model. It leads to asym- 
metric lateral domains due to unsynchronized winding 
of DW in the depth, and may clarify the present confu- 
sion and debate based on lateral multi-domain observa- 
tions. It also explains the long tail of the hysteresis loops, 
which would otherwise disappear if a$ is not included as 
in Kiwi's model (Fig. |3| inset) [28j [. Since a square hys- 
teresis loop is observed above T/v, this low temperature 
behavior must arise from the interfacial inhomogeneity. 
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FIG. 3: (color online) Simulation of the longitudinal (open 
symbols) and transverse (filled symbols) hysteresis loop con- 
sidering 50 kJ/m 3 in-plane uniaxial anisotropy and 0.5° 
fanning of AF uncompensated moment orientation. (Inset) 
Simulated transverse hysteresis loops with uniform (red open 
symbols) and inhomogeneous (black line) interfacial coupling. 

So far we demonstrated that the local IDWs nucleated 
in approaching the negative saturation cause the asym- 
metric reversal. This result is unambiguously confirmed 
by MOKE experiments probing the FM-air and FM-AF 
interfaces independently. In this experiment, a sample 
with MgF 2 (110) / FeF 2 (50 nm) / Py (70 nm) / Al 
(4 nm) is cooled below Tn in fio Hpc = 0.02 T, and 
MOKE is performed on both the top and bottom sur- 
faces of the sample with HeNe laser (A = 632.8 nm) at 
45 degree incidence (Fig. 21 inset (c)). Probing the depth 
dependence of the FM domain structure is possible be- 
cause the 28 nm penetration depth of the light |29[ is less 
than half of the Py thickness, and both MgF 2 and FeF 2 
are transparent. A clear difference is seen between the 
two MOKE measurements (Fig. Probing the FM- 
AF interface shows a much more rounded and longer tail 
compared with the one from FM-air interface, confirm- 
ing the existence of domain structures in the depth. The 
sample was also measured using SQUID magnetometry 
to which the entire sample contributes equally. The re- 
sultant hysteresis loop lies between the two MOKE loops. 

We also performed micromagnetic simulations under 
identical assumptions using the same parameters as 
above to generate the random frozen AF moments [3(j. 
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FIG. 4: (color online) Experiment (open symbol) and micro- 
magnetic simulation (solid line) on FeF2 (70 nm) / Py (70 nm) 
at 10 K after field cooling in a 0.02 T field. Experimental 
curves obtained from MOKE measurement from the FM-air 
(blue triangle) and FM-AF (red circle) interface and SQUID 
magnetometry (black square). The schematic of the MOKE 
experiment is shown in inset (c). The upper- left inset shows 
the simulated FM spin configuration (500x500 nm 2 ) at the 
FM-AF (a) and FM-air (b) interface at fi H = -0.06 T. 

The exponential decay of MOKE in the FM is simulated 
by giving each FM discretization layer in the depth an ap- 
propriate weight according to the Py 28 nm penetration 
depth. A very good agreement is obtained for all three 
hysteresis loops simultaneously with a slight adjustment 
of the inter facial coupling j^J. At n^H = -0.06 T, a large 
difference between the two MOKE measurements is ob- 
served. The simulated spin configuration at this field 
shows that the FM close to the FM-AF interface is only 
partially reversed forming lateral domain patterns, while 
at the FM-air interface the FM is fully reversed (inset (a) 
and (b) of Figure EJ. This confirms that the local IDW 
model leads to asymmetrically rounded hysteresis loops. 

In summary, we found strongly asymmetric hysteresis 
loops in a simple model exchange bias system FeF2/(Ni, 
Py). By combining vector magnetometry, MOKE with 
micromagnetic simulation, we clearly showed that the 
asymmetric reversal directly results from the FM domain 
structure in the depth due to the broken symmetry at the 
interface. The hotly debated issue over the asymmetric 
reversal process over the past 5 years solely focused on 
lateral FM domains, and its origin was controversial until 
now. FM parallel domains were predicted j^,!^. How- 
ever, they were not confirmed experimentally. They were 
mostly ignored in microscopy studies p| and simulations 
generally assu ming the FM to be a single moment or 
one monolayer . This situation was mostly due to the 
weak coupling at the FM/AF interface, and limitations of 
different experimental and modeling techniques. Disper- 
sions in AF crystallinity and anisotropy also smear out 
manifestations of parallel domain walls. Our study of a 



simple EB model system, combining different experimen- 
tal and simulation techniques, unambiguously demon- 
strates the presence of such domains and their dominant 
role on the asymmetric reversal. 
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